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SPEECH PROCESSING SYSTEM 



The present invention relates to an apparatus for and 
method of speech processing- The invention particularly 
relates to a statistical processing of an input speech 
signal to derive parameters defining the speech 
production system which generated the speech. 

Most speech processing systems, such as speech 
recognition systems, speech encoding systems, speaker 
verification systems and the like all include a signal 
processing front end, which converts the speech waveform 
into a parametric representation. One of the most common 
types of parameters (or coefficients) which are generated 
to represent speech are linear prediction (LP) 
coefficients. To represent the entire input speech 
signal, the speech is divided into a number of time 
frames (typically having a duration of 10-30ms) and a set 
of LP coefficients is calculated to represent the speech 
within each time frame. This LP analysis assumes that 
the value of a current speech sample can be predicted 
from a linear weighted combination of the k most recent 
speech samples. Based on this model, the task of the LP 
analysis is to identify the value of the weightings (or 
coefficients) which minimises the mean squared error 



between the actual value of the current speech sample and 
the predicted value of the current speech sample. 

One of the problems with this LP analysis is that it 
performs the analysis of the speech within each frame in 
isolation from the speech within other frames. It also 
assumes that the same number of weightings or 
coefficients will represent the speech within each time 
frame. As a result, errors can be introduced making 
subsequent recognition processing or the like more prone 
to error. 

An aim of the present invention is to provide an 
alternative technique for pre-possessing speech prior to 
use in a speech recognition system, a speech encoding 
system, a speaker verification system or the like. 

Exemplary embodiments of the present invention will now 
be described with reference to the accompanying drawings 
in which: 

Figure 1 is a schematic view of a computer which may be 
progreunmed to operate in accordance with an embodiment of 
the present invention; 



10. 
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Figure 2 is a block diagram illustrating the principal 
components of a speech recognition system; 

Figure 3 is a block diagram representing a model employed 
by a statistical analysis unit which forms part of the 
speech recognition system shown in Figure 2; 

Figure 4 is a flow chart illustrating the processing 
steps performed by a model order selection unit forming 
part of the statistical analysis unit shown in Figure 2; 

Figure 5 is a flow chart illustrating the main processing 
steps employed by a Simulation Smoother which forms part 
of the statistical analysis unit shown in Figure 2; 

Figure 6 is a block diagram illustrating the main 
processing components of the statistical analysis unit 
shown in Figure 2; 



20 Figure 7 is a memory map illustrating the data that is 

stored in a memory which forms part of the statistical 
analysis unit shown in Figure 2 ; 
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Figure 8 is a flow chart illustrating the main processing 
steps performed by the statistical analysis unit shown in 




4 

Figure 6; 

Figure 9a is a histogram for a model order of an auto 
regressive filter model which forms part of the model 
5 shown in Figure 3; 

Figure 9b is a histogram for the variance of process 
noise modelled by the model shown in Figure 3; 

10 Figure 9c is a histogram for a third coefficient of the 

AR filter model; 

Figure 10 is a block diagreim illustrating the main 
components of an alternative speech recognition system in 
15 which data output by the statistical analysis unit is 

used to detect the beginning and end of speech within the 
input signal; 

Figure 11 is a schematic block diagram illustrating the 
20 principal components of a speaker verification system; 

Figure 12 is a schematic block diagram illustrating the 
principal components of an acoustic classification 
system: 



Figure 13 is a schematic block diagram illustrating the 
principal components of a speech encoding and 
transmission, and 

Figure 14 is a block diagreim illustrating the principal 
components of a data file annotation system which uses 
the statistical analysis unit shown in Figure 6 to 
provide quality of speech data for an associated 
annotation . 

Embodiments of the present invention can be implemented 
on computer hardware, but the embodiment to be described 
is implemented in software which is run in conjunction 
with processing hardware such as a personal computer, 
workstation, photocopier, facsimile machine or the like. 

Figure 1 shows a personal computer (PC) 1 which may be 
programmed to operate an embodiment of the present 
invention. A keyboard 3, a pointing device 5, a 
microphone 7 and a telephone line 9 are connected to the 
PC 1 via an interface 11. The keyboard 3 and pointing 
device 5 allow the system to be controlled by a user. 
The microphone 7 converts the acoustic speech signal of 
the user into an equivalent electrical signal and 
supplies this to the PC 1 for processing. An internal 



modem and speech receiving circuit (not shown) may be 
connected to the telephone line 9 so that the PC 1 can 
communication with, for example, a remote computer or 
with a remote user. 

The progrcun instructions which make the PC 1 operate in 
accordance with the present invention may be supplied for 
use with an existing PC 1 on, for example, a storage 
device such as a magnetic disc 13, or by downloading the 
software from the Internet (not shown) via the internal 
modem and telephone line 9- 

The operation of a speech recognition system embodying 
the present invention will now be described with 
reference to Figure 2. Electrical signals representative 
of the input speech from the microphone 7 are input to a 
filter 15 which removes unwanted frequencies (in this 
embodiment frequencies above 8 kHz) within the input 
signal. The filtered signal is then sampled (at a rate 
of 16 kHz) and digitised by the analogue to digital 
converter 17 and the digitised speech samples are then 
stored in a buffer 19. Sequential blocks (or frames) of 
speech samples are then passed from the buffer 19 to a 
statistical analysis unit 21 which performs a statistical 
analysis of each frame of speech samples in sequence to 



determine, amongst other things, a set of auto regressive 
(AR) coefficients representative of the speech within the 
frame. In this embodiment, the AR coefficients output by 
the statistical analysis unit 21 are then input, via a 
coefficient converter 23 to a cepstral based speech 
recognition unit 25. In this embodiment, therefore, the 
coefficient converter 23 converts the AR coefficients 
output by the analysis unit 21 into cepstral 
coefficients. This can be achieved using the conversion 
technique described in, for example, "Fundcimentals of 
Speech Recognition" by Rabiner and Juang at pages 115 and 
116. The speech recognition unit 2 5 then compares the 
cepstral coefficients for successive frames of speech 
with a set of stored speech models 21, which may be 
template based or Hidden Markov Model based, to generate 
a recognition result. 

Statistical Analysis Unit - Theory and Overview 

As mentioned above, the statistical analysis unit 21 
analyses the speech within successive frames of the input 
speech signal. In most speech processing systems, the 
frames are overlapping. However, in this embodiment, the 
frames of speech are non-overlapping and have a duration 
of 20ms which, with the 16kHz sampling rate of the 
analogue to digital converter 17, results in a frame size 



8 



of 320 samples. 

In order to perform the statistical analysis on each of 
the frames, the analysis unit 21 assumes that there is an 
underlying process which generated each sample within the 
frame. The model of this process used in this embodiment 
is shown in Figure 3. As shown, the process is modelled 
by a speech source 31 which generates, at time t = n, a 
raw speech sample s(n)- Since there are physical 
constraints on the movement of the speech articulators, 
there is some correlation between neighbouring speech 
samples. Therefore, in this embodiment, the speech 
source 31 is modelled by an auto regressive (AR) process. 
In other words, the statistical analysis unit 21 assumes 
that a current raw speech sample (s(n)) can be determined 
from a linear weighted combination of the most recent 
previous raw speech samples, i.e.: 



s(n) = a^s(n-l) + a^sin-l) + + a^s(n-k) + e(n) (1) 

where aj, aa a^ are the AR filter coefficients 

representing the amount of correlation between the speech 
samples; k is the AR filter model order; and e(n) 
represents random process noise which is involved in the 
generation of the raw speech samples. As those skilled in 



the art of speech processing will appreciate, these AR 
filter coefficients are the same coefficients that the LP 
analysis estimates albeit using a different processing 
technique. 

As shown in Figure 3, the raw speech samples s(n) 
generated by the speech source are input to a channel 33 
which models the acoustic environment between the speech 
source 31 and the output of the analogue to digital 
converter 17. Ideally, the channel 33 should simply 
attenuate the speech as it travels from the source 31 to 
the microphone. However, due to reverberation and other 
distortive effects, the signal {y(n)) output by the 
analogue to digital converter 17 will depend not only on 
the current raw speech sample (s(n)) but it will also 
depend upon previous raw speech samples. Therefore, in 
this embodiment, the statistical analysis unit 21 models 
the channel 33 by a moving average (MA) filter, i.e.: 

y(ri) = /io5(w) + h^s(p-l) + h^s{n-2) + + h^s(n-r) + e(«) (2) 

where y(n) represents the signal sample output by the 
analogue to digital converter 17 at time t = n; hp, h^, 

hj hr are the channel filter coefficients representing 

the amount of reverberation within the channel 33; r is 
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the channel filter model order; and e(n) represents a 
random additive measurement noise component. 

For the current frsune of speech being processed, the 
filter coefficients for both the speech source and the 
channel are assumed to be constant but unknown - 
Therefore, considering all N samples (where N = 360) in 
the current frame being processed gives: 

s(n) = a^s(ri-l) + a^s(n-2) + + af^s(n~k) + e{n) 

s{n-\) = a^s{n-2) + a^s{n-3) + + a^s{n-k-\) + e(«-l) 



(3) 



s{n-N-^\) = a^s{n-N) + a^s{n-N-\) + + aj^s{n-k-N+\) + e{n-N+\) 



which can be written in vector form as 



= S.a + 



(4) 



20 



25 



where 



S = 



s{n-l) sin~2) sin-3) 
sin -2) sin -3) 5(«-4) 
sin-3) sin -4) 5(«-5) 



sin-N) sin-N-l) sin-N-2) . 



sin-k) 
sin-k-\) 
sin-k-2) 



sin-k-N+\) 



Nxk 



and 
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10 



i5 



20 



a = 



krl 



Sin) 
s{n-2) 



e{n-\) 
e{n-2) 



As will be apparent from the following discussion, it is 
also convenient to rewrite equation (3) in terms of the 
random error component (often referred to as the 
residual) e(n)- This gives: 

e(n) = s{n) - a^sin-l) - a^s(n'2) - - aj^s(n-k) 

e(w-l) = s(n-\) - a^sin-2) - a^s{n-3) - 

(5) 

e(n-N+\) = s(r2-N+\) - a^sin-N) - a^s{n-N-V) - - a^s{n-k-N+l) 

which can be written in vector notation as: 

= A^in) ( 6 ) 

where 



A = 



25 



1 


-^1 




-03 . . 


- 


0 


0 0 . 


. 0 


0 


1 






- -^k-\ 




0 0 . 


. 0 


0 


0 


1 


-a, . . 






-a^ 0 . 


. . 0 
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Similarly, considering the channel model defined by 
equation (2), with ho = 1 (since this provides a more 
stable solution), gives: 

q(n)=h^sin-\) +h2s{n-2) + +^^j(«-r) +e(«) 
5 qin-l) =h^s(ri-2) +h^sin'3)+.. +/i^5(n-r-l) +e(M-l) 

(7) 

q{n-N+\) =h^s{n-N) +h2s{n-N-\) ^ +h^s{n-r-N^l) ^z{n-N+\) 

(where q(n) = y(n) - s(n)) which can be written in vector 
10 form as: 

^(n) = Y,h (8) 

where 



15 



Y = 



^(«-l) s(n-2) sin -3) 
s{n-2) 5(n-3) j(«-4) 
s(n-3) s(n-4) sin -5) 



sin-r) 
sin-r-\) 
sin-r -2) 



sin-N) sin-N-l) sin'N-2) 



sin-r-N+l) 



Nxr 



20 



and 



h = 



rxl 



qin) 

qin-\) 
q(n-2) 



qin-N+l) 



A6f/ 



£(") = 



Bin) 
e(n-l) 

E(«-2) 



tin-N+\) 
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In this embodiment, the analysis unit 21 aims to 
determine, amongst other things, values for the AR filter 
coefficients (a) which best represent the observed signal 
samples (y(n)) in the current frame. It does this by 
determining the AR filter coefficients (a) that maximise 
the joint probability density function of the speech 
model, channel model, speech samples and the noise 
statistics given the observed signal samples output from 
the analogue to digital converter 17, i.e. by 
determining : 

{ p(.a,k.h.r.cl,oU(.n)\m) } 

where a^^ and a^^ represent the process and measurement 
noise statistics respectively- As those skilled in the 
art will appreciate, this function defines the 
probability that a particular speech model, channel 
model, raw speech samples and noise statistics generated 
the observed frame of speech samples (y:(n) ) from the 
analogue to digital converter- To do this, the 
statistical analysis unit 21 must determine what this 
function looks like. This problem can be simplified by 
rearranging this probability density function using Bayes 
law to give: 



/?Oi(^)U(A»),^,r,a:)/7(^(^)[fl,j^,q:);?(a|A)/?(^|r)/7(aJ/7(a,);?(/:)/?(r) ^^^^ 

POi(n)) 

As those skilled in the art will appreciate, the 
denominator of equation (10) can be ignored since the 
probability of the signals from the analogue to digital 
converter is constant for all choices of model. 
Therefore, the AR filter coefficients that maximise the 
function defined by equation (9) will also maximise the 
numerator of equation (10). 

Each of the terms on the numerator of equation (10) will 
now be considered in turn- 

P(s(n) I a, kr cr/; 

This term represents the joint probability density 
function for generating the vector of raw speech samples 
(s(n)) during a frame, given the AR filter coefficients 
(a), the AR filter model order (k) and the process noise 
statistics (Oe^). From equation (6) above, this joint 
probability density function for the raw speech samples 
can be determined from the joint probability density 
function for the process noise. In particular p(s(n) |a, 
k, Oe^) is given by: 



15 



(11) 



10 



where p<e(n)) is the joint probability density function 
for the process noise during a frame of the input speech 
and the second term on the right-hand side is known as 
the Jacobean of the transformation. In this case, the 
Jacobean is unity because of the triangular form of the 
matrix A ( see equations ( 6 ) above ) • 



20 



In this embodiment, the statistical analysis unit 21 
assumes that the process noise associated with the speech 
source 31 is Gaussian having zero mean and some unknown 
variance a^^. The statistical analysis unit 21 also 
assumes that the process noise at one time point is 
independent of the process noise at another time point. 
Therefore, the joint probability density function for the 
process noise during a frame of the input speech (which 
defines the probability of any given vector of process 
noise e{n) occurring) is given by: 



P(m) = (27101) ' exp 



2a 



(12) 



25 



Therefore, the joint probability density function for a 



vector of raw speech samples given the AR filter 
coefficients (a), the AR filter model order (k) and the 
process noise variance (Oq^) is given by: 



N 



p(s.(n)\a,k,aj = (27caJ ^ exp 



2a 



(13) 



p(y:(n)\s(n), r, a/) 

This term represents the joint probability density 
function for generating the vector of speech saonples 
(^(n) ) output from the analogue to digital converter 17, 
given the vector of raw speech samples (^(n) ) , the 
channel filter coefficients (h) , the channel filter model 
order (r) and the measurement noise statistics (a^^). 
From equation (8), this joint probability density 
function can be determined from the joint probability 
density function for the process noise- In particular, 
P(y(n)|s(n), h, r, a,^) is given by: 

where p(e.(n) ) is the joint probability density function 
for the measurement noise during a freune of the input 
speech and the second term on the right hand side is the 
Jacobean of the transformation which again has a value of 



6£(«) 



£(«) = ^(n) - Yh 



(14) 
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In this embodiment, the statistical analysis unit 21 
assumes that the measurement noise is Gaussian having 
zero mean and some unknown variance a^^ - It also assumes 
that the measurement noise at one time point is 
independent of the measurement noise at another time 
point. Therefore, the joint probability density function 
for the measurement noise in a frame of the input speech 
will have the same form as the process noise defined in 
equation (12). Therefore, the joint probability density 
function for a vector of speech samples (y(n) ) output 
from the analogue to digital converter 17, given the 
channel filter coefficients (h) , the channel filter model 
order (r), the measurement noise statistics {o^^) and the 
raw speech samples (s(n) ) will have the following form: 



2a 



(15) 



As those skilled in the art will appreciate, although 
this joint probability density function for the vector of 
speech samples (y:(n)) is in terms of the variable a(n), 
this does not matter since g(n) is a function of y(n) and 
s(n), and s(n) is a given variable (ie known) for this 
probability density function. 
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This term defines the prior probability density function 
for the AR filter coefficients (a) and it allows the 
5 statistical analysis unit 21 to introduce knowledge about 

what values it expects these coefficients will take. In 
this embodiment, the statistical analysis unit 21 models 
this prior probability density function by a Gaussian 
having an unknown variance (Oa^) and mean vector iUa) r 
10 i.e. : 



N 

p(a\k,a^^^ ) = (2kcI) ^ exp 



2ol 



(16) 



By introducing the new variables and the prior 

15 density functions {pio^^) and pilu)) for these variables 

must be added to the numerator of equation (10) above. 
Initially, for the first frame of speech being processed 
the mean vector (jju) can be set to zero and for the 
second and subsequent frames of speech being processed, 
20 it can be set to the mean vector obtained during the 

processing of the previous freune. In this case, p(iu) is 
just a Dirac delta function located at the current value 
of lu and can therefore be ignored. 
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With regard to the prior probability density function for 



the variance of the AR filter coefficients, the 
statistical analysis unit 21 could set this equal to some 
constant to imply that all variances are equally 
probable. However, this term can be used to introduce 
knowledge about what the variance of the AR filter 
coefficients is expected to be. In this embodiment, 
since variances are always positive, the statistical 
analysis unit 21 models this variance prior probability 
density function by an Inverse Gamma function having 
parameters and Pa, i.e.: 

At the beginning of the speech being processed, the 
statistical analysis unit 21 will not have much knowledge 
about the variance of the AR filter coefficients. 
Therefore, initially, the statistical analysis unit 21 
sets the variance Oa^ and the a and 3 parameters of the 
Inverse Gamma function to ensure that this probability 
density function is fairly flat and therefore non- 
informative. However, after the first frame of speech 
has been processed, these parameters can be set more 
accurately during the processing of the next frame of 
speech by using the parameter values calculated during 
the processing of the previous frame of speech. 



-1 



(17) 
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P(h\r) 

This term represents the prior probability density 
function for the channel model coefficients (h) and it 
allows the statistical analysis unit 21 to introduce 
5 knowledge about what values it expects these coefficients 

to take- As with the prior probability density function 
for the AR filter coefficients, in this embodiment, this 
probability density function is modelled by a Gaussian 
having an unknown variance (Oh^) and mean vector (Uh) , 
10 i.e.: 



_N 

pit\r,al^^) = (27Lal) ^ exp 



2al 



(18) 



15 Again, by introducing these new variables, the prior 

density functions (p(ah) and p((Mh)) must be added to the 
numerator of equation (10). Again, the mean vector can 
initially be set to zero and after the first frame of 
speech has been processed and for all subsequent frames 

20 of speech being processed, the mean vector can be set to 

equal the mean vector obtained during the processing of 
the previous frame. Therefore, p(Mh) is also just a 
Dirac delta function located at the current value of ^ 
and can be ignored. 
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With regard to the prior probability density function for 
the variance of the channel filter coefficients, again, 
in this embodiment, this is modelled by an Inverse Gamma 
function having parameters and Ph- Again, the variance 
5 (Oh^) and the a and (5 parameters of the Inverse Gamma 

function can be chosen initially so that these densities 
are non-informative so that they will have little effect 
on the subsequent processing of the initial frame - 

10 P(^J) and p(cr/) 

These terms are the prior probability density functions 
for the process and measurement noise variances and 
again, these allow the statistical analysis unit 21 to 
introduce knowledge about what values it expects these 

15 noise variances will take. As with the other variances, 

in this embodiment, the statistical analysis unit 21 
models these by an Inverse Gamma function having 
'parameters 3e and a^, 3, respectively. Again, these 

variances and these Gamma function parameters can be set 

20 initially so that they are non-informative and will not 

appreciably affect the subsequent calculations for the 
initial ^ frame . 
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p(k) and p(r) 

These terms are the prior probability density functions 




22 

for the AR filter model order (k) and the channel model 
order (r) respectively. In this embodiment, these are 
modelled by a uniform distribution up to some maximum 
order- In this way, there is no prior bias on the number 
5 of coefficients in the models except that they can not 

exceed these predefined maximums- In this embodiment, 
the maximum AR filter model order (k) is thirty and the 
maximum channel model order (r) is one hundred and fifty. 



10 



Therefore, inserting the relevant equations into the 
numerator of equation (10) gives the following joint 
probability density function which is proportional to 
P(a,k,h,r,aa^ah2,ae2,a,2,s(n) |y(n) ) : 



15 



20 



N 

2^2, 



(27cag) 2 exp 



-I 
2a? 



N 

.2^ 2, 



X (2nOg) exp 



-1 



lot 



_N 

X (2nol) 2 exp 



■(fl-M/(a-H; 



2o 



N 

2\ 2 



X (InOfj) 2 exp 



2af 



X exp 


-1 


X exp 

P*r(a,) 


-1 


.°«p<.. 




X exp 

Per(a.) 


-1 

g^p. 


X exp 

Per(a,) 


-1 



(19) 
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Gibbs Sampler 

In order to determine the form of this joint probability 
density function, the statistical analysis unit 21 "draws 
samples" from it- In this embodiment, since the joint 
5 probability density function to be sampled is a complex 

multivariate function, a Gibbs sampler is used which 
breaks down the problem into one of drawing ssunples from 
probability density functions of smaller dimensionality. 
In particular, the Gibbs sampler proceeds by drawing 
10 random variates from conditional densities as follows: 

first iteration 

p(a,^/2°,r^af,af,af,af,£(/i)VW) ^ a\k' 

p(of|a^^^Zl^r^af,af,0J',;£(/l)V«)) ot 

second iteration 
^° p(aA16',r'.o^',af,oj'.£(/»)'j^/j)) -* ^,k^ 

/7(Ma'.*'.cj^'.''c'.cJa.cJ*'X«)'j'(")) AV" 
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etc . 
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where ( h% r\ {o,^)\ ioJ)\ {o^^)\ s(n)°) are 

initial values which may be obtained from the results of 
the statistical analysis of the previous freune of speech, 
or where there are no previous frcunes, can be set to 
5 appropriate values that will be known to those skilled in 

the art of speech processing. 
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As those skilled in the art will appreciate, these 
conditional densities are obtained by inserting the 
current values for the given (or known) variables into 
the terms of the density function of equation (19). For 
the conditional density p(a,k|...) this results in: 



p(a,k\...) « exp 



— {&inY&(n) - 2a^Ss.(n) + a^S'^S a) 



2a! 



X exp 



lot 



(20) 



which can be simplified to give: 
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p(a,Al..-) « exp 



2 



T 



+ a' 



+ — 

2 2 



/ J 



(21) 
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which is in the form of a standard Gaussian distribution 
having the following covariance matrix: 
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2 2 



-I 



(22) 



The mean value of this Gaussian distribution can be 
determined by differentiating the exponent of equation 
(21) with respect to a and determining the value of a 
which makes the differential of the exponent equal to 
zero- This yields a mean value of: 





-1 




2 2 




2 



(23) 
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A sample can then be drawn from this standard Gaussian 
distribution to give a^ (where g is the g^^ iteration of 
the Gibbs sampler) with the model order (k^) being 
determined by a model order selection routine which will 
be described later- The drawing of a sample from this 
Gaussian distribution may be done by using a random 
number generator which generates a vector of random 
values which are uniformly distributed and then using a 
transformation of random variables using the covariance 
matrix and the mean value given in equations (22) and 
(23) to generate the sample. In this embodiment, 
however, a random number generator is used which 
generates random numbers from a Gaussian distribution 
having zero mean and a variance of one- This simplifies 
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the transformation process to one of a simple scaling 
using the covariance matrix given in equation (22) and 
shifting using the mean value given in equation (23). 
Since the techniques for drawing samples from Gaussian 
distributions are well known in the art of statistical 
analysis, a further description of them will not be given 
here. A more detailed description and explanation can be 
found in the book entitled "Numerical Recipes in C", by 
W. Press et al, Cambridge University Press, 1992 and in 
particular at chapter 7 . 
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As those skilled in the art will appreciate, however, 
before a scimple can be drawn from this Gaussian 
distribution, estimates of the raw speech Scimples must be 
available so that the matrix S and the vector s(n) are 
known- The way in which these estimates of the raw 
speech samples are obtained in this embodiment will be 
described later - 
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A similar analysis for the conditional density p(h,r| - . . ) 
reveals that it also is a standard Gaussian distribution 
but having a covariance matrix and mean value given by: 
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(24) 
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from which a sample for h^ can be drawn in the manner 
described above, with the channel model order (r^) being 
determined using the model order selection routine which 
will be described later. 

5 

A similar analysis for the conditional density p(ae^|---) 
shows that: 

Pioll- ) - (cT')'"exp 

10 

where : 

E = ^(nf^in) - 2a^Ss(n) + a^S^S a 



-E 
2ol 



exp 



(25) 



which can be simplified to give: 



exp 



-1 



E 
2 



(26) 
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which is also an Inverse Gamma distribution having the 
following parameters: 
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(27) 
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A sample is then drawn from this Inverse Gamma 
distribution by firstly generating a random number from 
a uniform distribution and then performing a 
transformation of random variables using the alpha and 
beta parameters given in equation (27), to give (Oe^)^. 

A similar analysis for the conditional density p{o^^\ ...) 
reveals that it also is an Inverse Gamma distribution 
having the following parcuneters: 



~ ^ Pe = (28) 



where : 



15 E* = q.(ny<i{n) - 2h^ Yg.(ri) + b^Y^Vh 

A sample is then drawn from this Inverse Gamma 
distribution in the manner described above to give (ae^)^. 

20 A similar analysis for conditional density pCOa^l,..) 

reveals that it too is an Inverse Gamma distribution 
having the following parcimeters : 



N - 2P^ 

«a = T K = ; (29) 

25 - 



A sample is then drawn from this Inverse Gamma 
distribution in the manner described above to give (a^^)?. 



Similarly, the conditional density piOhM--*) is also an 
Inverse Gamma distribution but having the following 
parameters : 

d. = — + a . and = (30) 

' 2 ' '2. p,.(&-|i/(^-|i,) 

A sample is then drawn from this Inverse Gamma 
distribution in the manner described above to give (ah^)^. 



As those skilled in the art will appreciate, the Gibbs 
sampler requires an initial transient period to converge 
to equilibrium (known as burn-in). Eventually, after L 
iterations, the sample (a^, kS h^, r^, (Oe^)^, {o^^)^, 
(Oa^)^, (Oh^)^, s(n)^) is considered to be a sample from the 
joint probability density function defined in equation 
(19). In this embodiment, the Gibbs sampler performs 
approximately one hundred and fifty (150) iterations on 
each frcime of input speech and discards the Scimples from 
the first fifty iterations and uses the rest to give a 
picture (a set of histograms) of what the joint 
probability density function defined in equation (19) 
looks like- From these histogreims, the set of AR 
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coefficients (a) which best represents the observed 
speech samples (y(n) ) from the analogue to digital 
converter 17 are determined. The histograms are also 
used to determine appropriate values for the variances 
5 and channel model coefficients (h) which can be used as 

the initial values for the Gibbs seunpler when it 
processes the next frame of speech. 

Model Order Selection 

10 As mentioned above, during the Gibbs iterations, the 

model order (k) of the AR filter and the model order (r) 
of the channel filter are updated using a model order 
selection routine. In this embodiment, this is performed 
using a technique called "Reversible jump Markov chain 

15 Monte Carlo computation". This technique is described in 

the paper entitled "Reversible jump Markov chain Monte 
Carlo Computation and Bayesian model determination" by 
Peter Green, Biometrika, vol 82, pp 711 to 732, 1995. 

2 0 Figure 3 is a flow chart which illustrates the processing 

steps performed during this model order selection routine 
for the AR filter model order (k) . As shown, in step si, 
a new model order (kj) is proposed. In this embodiment, 
the new model order will normally be proposed as k2 = ki 

25 ± Ir but occasionally it will be proposed as k2 = ki ± 2 
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and very occasionally as = ki ± 3 etc. To achieve 
this, a sample is drawn from a discretised Laplacian 
density function centred on the current model order (kj) 
and with the variance of this Laplacian density function 
being chosen a priori in accordance with the degree of 
sampling of the model order space that is required. 

The processing then proceeds to step s3 where a model 
order variable (MO) is set equal to: 



MO = max < 



(31) 



where the ratio term is the ratio of the conditional 
probability given in equation (21) evaluated for the 

15 current AR filter coefficients (a) drawn by the Gibbs 

sampler for the current model order (ki) and for the 
proposed new model order (kj). If kj > ki, then the 
matrix S must first be resized and then a new sample 
must be drawn from the Gaussian distribution having the 

20 mean vector and covariance matrix defined by equations 

(22) and (23) (determined for the resized matrix S),to 
provide the AR filter coefficients (a<i,,c2>) for the new 
model order (ks). If kz < ki then all that is required is 
to delete the last (ki - kj) Scunples of the a vector. If 

2 5 the ratio in equation (31) is greater than one, then this 
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implies that the proposed model order (kj) is better than 
the current model order whereas if it is less than one 
then this implies that the current model order is better 
than the proposed model order. However, since 

5 occasionally this will not be the case, rather than 

deciding whether or not to accept the proposed model 
order by comparing the model order variable (MO) with a 
fixed threshold of one, in this embodiment, the model 
order variable (MO) is compared, in step s5, with a 

10 random number which lies between zero and one. If the 

model order variable (MO) is greater than this random 
number, then the processing proceeds to step s7 where the 
model order is set to the proposed model order (kz) and 
a count associated with the value of kz is incremented. 

15 If , on the other hand, the model order variable (MO) is 

smaller than the random number, then the processing 
proceeds to step s9 where the current model order is 
maintained and a count associated with the value of the 
current model order (ki) is incremented. The processing 

20 then ends. 

This model order selection routine is carried out for 
both the model order of the AR filter model and for the 
model order of the channel filter model. This routine 
2 5 may be carried out at each Gibbs iteration. However, 



this is no-t essential- Therefore, in this embodiment, 
this model order updating routine is only carried out 
every third Gibbs iteration. 

Simulatlan Smoother 

As mentioned above, in order to be able to draw samples 
using the Gibbs sampler, estimates of the raw speech 
samples are required to generate s{n), S and Y which are 
used in the Gibbs calculations. These could be obtained 
from the conditional probability density function 
p(s(n)|..0* However, this is not done in this 
embodiment because of the high dimensionality of S(n) . 
Therefore, in this embodiment, a different technique is 
used to provide the necessary estimates of the raw speech 
samples. In particular, in this embodiment, a 

"Simulation Smoother" is used to provide these estimates. 
This Simulation Smoother was proposed by Piet de Jong in 
the paper entitled "The Simulation Smoother for Time 
Series Models", Biometrika (1995), vol 82,2, pages 339 to 
350. As those skilled in the art will appreciate, the 
Simulation Smoother is run before the Gibbs Sampler. It 
is also run again during the Gibbs iterations in order to 
update the estimates of the raw speech samples. In this 
embodiment, the Simulation Smoother is run every fourth 
Gibbs iteration- 
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In order to run the Simulation Smoother, the model 
equations defined above in equations (4) and (6) must be 
written in "state space" format as follows: 



f(n) =/.£(?!- 1) + 



(32) 
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where 



A = 



flj flj ^3 

1 0 0 
0 1 0 



a, 0 

0 0 
0 0 



1 0 
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and 



sin) 
Hn-1) 
sin -2) 



sin-r+l) 



rxl 



tin) = 



ein) 
0 
0 



J rxl 



20 



25 



With this state space representation, the dimensionality 
of the raw speech vectors (£(n) ) and the process noise 
vectors (e(n) ) do not need to be Nxl but only have to be 
as large as the greater of the model orders - k and r. 
Typically, the channel model order (r) will be larger 
than the AR filter model order (k) - Hence, the vector of 
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raw speech samples (s(n)) and the vector of process noise 
(e(n)) only need to be rxl and hence the dimensionality 
of the matrix A only needs to be rxr. 

5 The Simulation Smoother involves two stages - a first 

stage in which a Kalman filter is run on the speech 
samples in the current frame and then a second stage in 
' which a "smoothing" filter is run on the speech samples 

in the current frame using data obtained from the Kalman 

10 filter stage. Figure 4 is a flow chart illustrating the 

processing steps performed by the Simulation Smoother. 
As shown, in step s21, the system initialises a time 
variable t to equal one. During the Kalman filter stage, 
this time variable is run from t = 1 to N in order to 

15 process the N speech samples in the current frame being 

processed in time sequential order. After step s21, the 
processing then proceeds to step s23, where the following 
Kalman filter equations are computed for the current 
speech sample (y(t)) being processed: 

kit) = (APm)-^ifr' 

^ (33) 

f(/+l) = Ai(t) k it) MO 



25 



L(0 - A - k^{t).h' 
+ = AP{t)L{tf + a^/ 
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where the initial vector of raw speech samples (§(1)) 
includes raw speech samples obtained from the processing 
of the previous frame (or if there are no previous frames 
then s(i) is set equal to zero for i < 1); P(l) is the 
5 variance of s(l) (which can be obtained from the previous 

frame or initially can be set to a^^); h is the current 
set of channel model coefficients which can be obtained 
from the processing of the previous freime (or if there 
are no previous frames then the elements of h can be set 

10 to their expected values - zero); y(t) is the current 

speech sample of the current frame being processed and I 
is the identity matrix. The processing then proceeds to 
step s25 where the scalar values w(t) and d(t) are stored 
together with the rxr matrix L(t) (or alternatively the 

15 Kalman filter gain vector k£(t) could be stored from 

which L(t) can be generated). The processing then 
proceeds to step s2 7 where the system determines whether 
or not all the speech samples in the current frame have 
been processed- Xf they have not^ then the processing 

20 proceeds to step s29 where the time variable t is 

incremented by one so that the next sample in the current 
frame will be processed in the same way. Once all N 
samples in the current frame have been processed in this 
way and the corresponding values stored, the first stage 

25 of the Simulation Smoother is complete. 



The processing then proceeds to step s31 where the second 
stage of the Simulation Smoother is started in which the 
smoothing filter processes the speech samples in the 
current frame in reverse sequential order. As shown, in 
step s31 the system runs the following set of smoothing 
filter equations on the current speech sample being 
processed together with the stored Kalman filter 
variables computed for the current speech sample being 
processed: 

C(0 = a^{/ - olUiO) 

u(0 ~mc(0) 

(34) 

u{t-\) - hd{ty'h^ L{tfu{t)L{t) f v{tfc{ty'v(t) 

i{t) = cliit) + n(0 where = [c(0 c{f-2) ... e{t-r+l)f 

m = Ai{t'\) + m where Ht) = [s{t) s\t-\) f(r-2) ... f^-r+l)]^ 

and lit) = [e(t) 0 0 ... 0 ]^ 

where ii(t) is a sample drawn from a Gaussian distribution 
having zero mean and covariance matrix C(t); the initial 
vector r(t=N) and the initial matrix U(t=N) are both set 
to zero; and s{0) is obtained from the processing of the 
previous frame (or if there are no previous frames can be 
set equal to zero) . The processing then proceeds to step 
s3 3 where the estimate of the process noise (e(t)) for 



38 

the current speech sample being processed and the 
estimate of the raw speech sample (s(t)) for the current 
speech sample being processed are stored- The processing 
then proceeds to step s35 where the system determines 
whether or not all the speech seimples in the current 
frame have been processed- If they have not, then the 
processing proceeds to step s37 where the time variable 
t is decremented by one so that the previous sample in 
the current frame will be processed in the same way. 
Once all N samples in the current frame have been 
processed in this way and the corresponding process noise 
and raw speech samples have been stored, the second stage 
of the Simulation Smoother is complete and an estimate of 
s(n) will have been generated. 

As shown in equations ( 4 ) and ( 8 ) , the matrix S and the 
matrix Y require raw speech samples s{n-N-l) to s(n-N- 
k+1) and s(n-N-l) to s(n-N-r+l) respectively in addition 
to those in s.(n) . These additional raw speech samples 
can be obtained either from the processing of the 
previous freune of speech or if there are no previous 
frames, they can be set to zero- With these estimates of 
raw speech Scunples, the Gibbs seunpler can be run to draw 
samples from the above described probability density 
functions . 
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Statistical Analysis Unit - Operation 

A description has been given above of the theory 
underlying the statistical analysis unit 21. A 
description will now be given with reference to Figures 
5 6 to 8 of the operation of the statistical analysis unit 

21 that is used in the embodiment - 

) Figure 6 is a block diagram illustrating the principal 

components of the statistical analysis unit 21 of this 

10 embodiment. As shown, it comprises the above described 

Gibbs sampler 41, Simulation Smoother 4 3 (including the 
Kalman filter 43-1 and smoothing filter 43-2) and model 
order selector 45- It also comprises a memory 47 which 
receives the speech samples of the current freime to be 

15 processed, a data analysis unit 4 9 which processes the 

data generated by the Gibbs sampler 41 and the model 
order selector 45 and a controller 50 which controls the 
operation of the statistical analysis unit 21. 

20 As shown in Figure 6, the memory 47 includes a non 

volatile memory area 47-1 and a working memory area 47-2. 
The non volatile memory 4 7-1 is used to store the joint 
probability density function given in equation (19) above 
and the equations for the variances and mean values and 

25 the equations for the Inverse Gamma parameters given 
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above in equations (22) to (24) and (27) to (30) for the 
above mentioned conditional probability density functions 
for use by the Gibbs sampler 41- The non volatile memory 
47-1 also stores the Kalman filter equations given above 
5 in equation (33) and the smoothing filter equations given 

above in equation 34 for use by the Simulation Smoother 
43. 

Figure 7 is a schematic diagram illustrating the 

10 parcimeter values that are stored in the working memory 

area (RAM) 47-2. As shown, the RAM includes a store 51 
for storing the speech samples yf ( 1 ) to yf (N) output by 
the analogue to digital converter 17 for the current 
frame (f) being processed. As mentioned above, these 

15 speech samples are used in both the Gibbs sampler 41 and 

the Simulation Smoother 43. The RAM 47-2 also includes 
a store 53 for storing the initial estimates of the model 
parameters (g=0) and the M scumples (g = 1 to M) of each 
parameter drawn from the above described conditional 

20 probability density functions by the Gibbs Scimpler 41 for 

the current frame being processed. As mentioned above, 
in this embodiment, M is 20 0 since the Gibbs sampler 41 
performs 250 iterations on each frame of input speech 
with the first fifty samples being discarded. The RAM 

25 47-2 also includes a store 55 for storing W(t), d(t) and 



L(t) for t = 1 to N which are calculated during the 
processing of the speech samples in the current frcune of 
speech by the above described Kalman filter 43-1. The 
RAM 47-2 also includes a store 57 for storing the 
estimates of the raw speech samples (Sf(t)) and the 
estimates of the process noise (ef(t)) generated by the 
smoothing filter 43-2, as discussed above. The RAM 47-2 
also includes a store 59 for storing the model order 
counts which are generated by the model order selector 45 
when the model orders for the AR filter model and the 
channel model are updated. 

Figure 8 is a flow diagreim illustrating the control 
program used by the controller 50, in this embodiment, to 
control the processing operations of the statistical 
analysis unit 21. As shown, in step s41, the controller 
50 retrieves the next frame of speech scimples to be 
processed from the buffer 19 and stores them in the 
memory store 51. The processing then proceeds to step 
s43 where initial estimates for the channel model, raw 
speech samples and the process noise and measurement 
noise statistics are set and stored in the store 53. 
These initial estimates are either set to be the values 
obtained during the processing of the previous frame of 
speech or, where there are no previous frames of speech. 
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are set to their expected values (which may be zero) . 
The processing then proceeds to step s45 where the 
Simulation Smoother 43 is activated so as to provide an 
estimate of the raw speech samples in the manner 
5 described above. The processing then proceeds to step s4 7 

where one iteration of the Gibbs sampler 41 is run in 
order to update the channel model, speech model and the 
process and measurement noise statistics using the raw 
speech samples obtained in step s45- These updated 
10 pareuneter values are then stored in the memory store 53- 

The processing then proceeds to step s49 where the 
controller 50 determines whether or not to update the 
model orders of the AR filter model and the channel 

15 model. As mentioned above, in this embodiment, these 

model orders- are updated every third Gibbs iteration. If 
the model orders are to be updated, then the processing 
proceeds to step s51 where the model order selector 45 is 
used to update the model orders of the AR filter model 

2 0 and the channel model in the manner described above. If 

at step s49 the controller 5 0 determines that the model 
orders are not to be updated, then the processing skips 
step s51 and the processing proceeds to step s53. At 
step s53, the controller 50 determines whether or not to 

25 perform another Gibbs iteration- If another iteration is 



43 

to be performed, then the processing proceeds to decision 
block s55 where th^ controller 50 decides whether or not 
to update the estimates of the raw speech samples (s(t))- 
If the raw speech samples are not to be updated, then the 
5 processing returns to step s47 where the next Gibbs 

iteration is run. 

As mentioned above, in this embodiment, the Simulation 
Smoother 4 3 is run every fourth Gibbs iteration in order 

10 to update the raw speech samples. Therefore, if the 

controller 50 determines, in step s55 that there has been 
four Gibbs iterations since the last time the speech 
samples were updated, then the processing returns to step 
s45 where the Simulation Smoother is run again to provide 

15 new estimates of the raw speech samples (s(t)). Once the 

controller 50 has determined that the required 150 Gibbs 
iterations have been performed, the controller 50 causes 
the processing to proceed to step s57 where the data 
analysis unit 49 analyses the model order counts 

20 generated by the model order selector 45 to determine the 

model orders for the AR filter model and the channel 
model which best represents the current frame of speech 
being processed. The processing then proceeds to step 
s59 where the data analysis unit 49 analyses the samples 

25 drawn from the conditional densities by the Gibbs sampler 
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41 to determine the AR filter coefficients (a), the 
channel model coefficients (h) , the variances of these 
coefficients and the process and measurement noise 
variances which best represent the current freime of 
5 speech being processed. The processing then proceeds to 

step s61 where the controller 50 determines whether or 
not there is any further speech to be processed. If 
there is more speech to be processed, then processing 
returns to step S41 and the above process is repeated for 
10 the next frame of speech. Once all the speech has been 

processed in this way, the processing ends. 

Data Analysis unit 

A more detailed description of the data analysis unit 49 
15 will now be given with reference to Figure 9. As 

mentioned above, the data analysis unit 49 initially 
determines, in step s57, the model orders for both the AR 
filter model and the channel model which best represents 
the current frame of speech being processed. It does 
20 this using the counts that have been generated by the 

model order selector 45 when it was run in step s51. 
These counts are stored in the store 59 of the RAM 47-2. 
In this embodiment, in determining the best model orders, 
the data analysis unit 49 identifies the model order 
25 having the highest count. Figure 9a is an exemplary 
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histogram which illustrates the distribution of counts 
that is generated for the model order (k) of the AR 
filter model. Therefore, in this example, the data 
analysis unit 49 would set the best model order of the AR 
5 filter model as five. The data analysis unit 49 performs 

a similar analysis of the counts generated for the model 
order (r) of the channel model to determine the best 
model order for the channel model. 

10 Once the data analysis unit 49 has determined the best 

model orders (k and r) , it then analyses the samples 
generated by the Gibbs sampler 41 which are stored in the 
store 53 of the RAM 47-2, in order to determine parameter 
values that are most representative of those samples. 

15 It does this by determining a histogram for each of the 

parameters from which it determines the most 
representative parameter value. To generate the 

histogram, the data analysis unit 49 determines the 
maximum and minimum sample value which was drawn by the 

20 Gibbs sampler and then divides the range of parameter 

values between this minimum and maximum value into a 
predetermined number of sub-ranges or bins. The data 
analysis unit 49 then assigns each of the sample values 
into the appropriate bins and counts how many samples 

25 are allocated to each bin. It then uses these counts to 



calculate a weighted average of the samples (with the 
weighting used for each sample depending on the count for 
the corresponding bin), to determine the most 
representative parameter value (known as the minimum mean 
square estimate (MMSE) ) . Figure 9b illustrates an 

example histogram which is generated for the variance 
(Oe^) of the process noise, from which the data analysis 
unit 4 9 determines that the variance representative of 
the sample is 0.3149. 

In determining the AR filter coefficients (ai for i = i 
to k) , the data analysis unit 49 determines and analyses 
a histogram of the samples for each coefficient 
independently. Figure 9c shows an exemplary histogram 
obtained for the third AR filter coefficient (aa), from 
which the data analysis unit 49 determines that the 
coefficient representative of the samples is -0.4977. 

In this embodiment, the data analysis unit 49 only 
outputs the AR filter coefficients which are passed to 
the coefficient converter 2 3 shown in Figure 2. The 
remaining parameter values determined by the data 
analysis unit 49 are stored in the RAM 47-2 for use 
during the processing of the next frame of speech. As 
mentioned above, the AR filter coefficients output by the 
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statistical analysis unit 21 are input to the coefficient 
converter 23 which converts these coefficients into 
cepstral coefficients which are then compared with stored 
speech models 27 by the speech recognition unit 25 in 
5 order to generate a recognition result. 

As the skilled reader will appreciate, a speech 
processing technique has been described above which uses 
statistical analysis techniques to determine sets of AR 

10 filter coefficients representative of an input speech 

signal. The technique is more robust and accurate than 
prior art techniques which employ maximum likelihood 
estimators to determine the AR filter coefficients. This 
is because the statistical analysis of each frame uses 

15 knowledge obtained from the processing of the previous 

frame. In addition, with the analysis performed above, 
the model order for the AR filter model is not assumed to 
be constant and can vary from frame to frame. In this 
way, the optimum number of AR filter coefficients can be 

20 used to represent the speech within each frame. As a 

result, the AR filter coefficients output by the 
statistical analysis unit 21 will more accurately 
represent the corresponding input speech. Further still, 
since the underlying process model that is used separates 

2 5 the speech source from the channel, the AR filter 
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coefficients that are determined will be more 
representative of the actual speech and will be less 
likely to include distortive effects of the channel. 
Further still, since variance information is available 
for each of the parameters, this provides an indication 
of the confidence of each of the parameter estimates. 
This is in contrast to, for example, linear prediction 
analysis where point estimates of the parameter values 
are determined. 

Alternative Embodiments 

In the above embodiment, the statistical analysis unit 
was used as a pre-processor for a speech recognition 
system in order to generate AR coefficients 
representative of the input speech. It also generated a 
number of other parameter values (such as the process 
noise variances and the channel model coefficients), but 
these were not output by the statistical analysis unit. 
As those skilled in the art will appreciate, the AR 
coefficients and some of the other parameters which are 
calculated by the statistical analysis unit can be used 
for other purposes. For example. Figure 10 illustrates 
a speech recognition system which is similar to the 
speech recognition system shown in Figure 2 except that 
there is no coefficient converter since the speech 
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recognition unit 2 5 and speech models 2 7 are AR 
coefficient based. The speech recognition system shown 
in Figure 10 also has an additional speech detection unit 
61 which receives the AR filter coefficients (a) together 
5 with the AR filter model order (k) generated by the 

statistical analysis unit 21 and which is operable to 
determine from them when speech is present within the 
signal received from the microphone 7. It can do this, 
since the AR filter model order and the AR filter 

10 coefficient values will be larger during speech than when 

there is no speech present. Therefore, by comparing the 
AR filter model order (k) and/or the AR filter 
coefficient values with appropriate threshold values, the 
speech detection unit 61 can determine whether or not 

15 speech is present within the input signal. When the 

speech detection unit 61 detects the presence of speech, 
it outputs an appropriate control signal to the speech 
recognition unit 25 which causes it to start processing 
the AR coefficients it receives from the statistical 

20 analysis unit 21. Similarly, when the speech detection 

unit 61 detects the end of speech, it outputs an 
appropriate control signal to the speech recognition unit 
25 which causes it to stop processing the AR coefficients 
it receives from the statistical analysis unit 21. 



25 
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In the above embodiments, a speech recognition system was 
described having a particular speech pre-processing front 
end which performed a statistical analysis of the input 
speech. As the those skilled in the art will appreciate, 
5 this pre-processing can be used in speech processing 

systems other than speech recognition systems . For 
example, as shown in Figure 11, the statistical analysis 
unit 21 may form a front end to a speaker verification 
system 65. In this embodiment, the speaker verification 
10 system 65 compares the sequences of AR filter 

coefficients output by the statistical analysis unit 21 
with pre-stored speaker models 67 to determine whether or 
not the received speech corresponds to a known user. 

15 Figure 12 illustrates another application for the 

statistical analysis unit 21. In particular. Figure 12 
shows an acoustic classification system. The statistical 
analysis unit 21 is used to generate the AR filter 
coefficients for an input audio signal (which may or may 

20 not be speech) in the manner described above. These 

coefficients are then passed to an acoustic 
classification system 66 which compares the AR 
coefficients from the statistical analysis unit with pre- 
stored acoustic models 6 8 to generate the classification 

25 result. This classification system may be used, for 
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example, to identify and categorise various input sounds, 
for example, percussion sounds, woodwind sounds, brass 
sounds as well as speech. 

5 Figure 13 illustrates another application for the 

statistical analysis unit 21. In particular. Figure 13 
shows a speech encoding and transmission system. The 
statistical analysis unit 21 is used to generate the AR 
filter coefficients in the manner described above, which 

10 are then passed to a channel encoder which encodes the 

sequences of AR filter coefficients so that they are in 
a more suitable form for transmission through a 
communications channel. The encoded AR filter 

coefficients are then passed to a transmitter 73 where 

15 the encoded data is used to modulate a carrier signal 

which is then transmitted to a remote receiver 75. The 
receiver 75 demodulates the received signal to recover 
the encoded data which is then decoded by a decoder 75. 
The sequences of AR filter coefficients output by the 

2 0 decoder are then either passed to a speech recognition 

unit 77 which compares the sequences of AR filter 
coefficients with stored reference models (not shown) to 
generate a recognition result or to a speech synthesis 
unit 79 which re-generates the speech and outputs it via 

2 5 a loud speaker 81. As shown, prior to application to the 



speech synthesis unit 79, the sequences, of AR filter 
coefficients may also pass through an optional processing 
unit 83 (shown in phantom) which can be used to 
manipulate the characteristics of the speech that is 
synthesised. One of the significant advantages of using 
the statistical analysis unit described above is that the 
model order for the AR filter model is not assumed to be 
constant and will vary from frame to frame. In this way, 
the optimum number of AR filter coefficients will be used 
to represent the speech within each frame. In contrast, 
with linear prediction analysis, the number of AR filter 
coefficients is assumed to be constant and hence the 
prior art techniques tend to over parameterise the speech 
in order to ensure that information is not lost. As a 
result, with the statistical analysis described above, 
the amount of data which has to be transmitted from the 
transmitter to the receiver will be less than with the 
prior art systems which assume a fixed size of AR filter 
model . 

Figure 14 shows another system which uses the statistical 
analysis unit 21 described above. The system shown in 
Figure 14. allows a user to add a voice annotation to a 
data file 91 for use in subsequent voice retrieval 
operations. In use, the user selects a data file to be 



annotated (which may be a photograph or the like). The 
user then speaks the voice annotation and the input 
speech received via the microphone 7 is processed by the 
statistical analysis unit 21 in the same way as described 
above, except the data analysis unit only outputs the 
variance {o^^) of the AR filter coefficients for each 
frcone of speech which is passed to a speech quality 
assessor 93. The speech quality assessor 93 uses this 
variance information as an indication of the quality of 
the input speech (although it will also depend on other 
factors). However, in general, when the input speech is 
of a high quality (i.e. not disturbed by high levels of 
background noise), the variance should be small and when 
there are high levels of noise, the variance should be 
large. The speech quality assessor 93 then outputs to a 
data file annotation unit 95 a signal identifying whether 
or not the input voice annotation is of a high quality. 
In this embodiment, the input speech is also passed 
through a traditional speech pre-processor 95 and then to 
a traditional speech recognition unit 97 which generates 
words and/or phonemes representative of the input speech 
which are then passed to the data file annotation unit 
99. The data file annotation unit 99 then annotates the 
selected data file 91 with the data output by the speech 
recognition unit 9 7 and the data output by the speech 
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quality assessor 93 and- stores the annotated data file in 
a database 101. 

As the those skilled in the art will appreciate, this 
5 speech quality indication which is stored with the data 

file is useful for subsequent retrieval operations. In 
particular, when the user wishes to retrieve a data file 
91 from the database 101 (using a voice query), it is 
useful to know the quality of the speech that was used to 

10 annotate the data file and/or the quality of the voice 

retrieval query used to retrieve the data file, since 
this will affect the retrieval perf onnance . In 
particular if the voice annotation is of a high quality 
and the user's retrieval query is also of a high quality, 

15 then a stringent search of the database 101 can be 

performed, in order to reduce the amount* of false 
identifications. In contrast, if the original voice 
annotation is of a low quality or if the user's retrieval 
query is of a low quality, then a less stringent search 

20 of the database 101 can be performed to give a higher 

chance of retrieving the correct data file 91. 



25 



In addition to using the variance of the AR filter 
coefficients as an indication of the speech quality, the 
variance (o^^) of the process noise is also a good measure 



of the quality of the input speech, since this variance 
is also measure of the energy in the process noise. 
Therefore, the variance of the process noise can be used 
in addition to or instead of the variance of the AR 
filter coefficients to provide the measure of quality of 
the input speech- 

In the above embodiment, a separate preprocessor 95 was 
provided as the front end to the speech recognition unit 
97. As those skilled in the art will appreciate, the 
statistical analysis unit 21 may act as the preprocessor 
for the speech recognition unit 97. Likewise, in the 
speech recognition system described above with reference 
to Figure 10, the statistical analysis unit 21 may be 
used solely for providing information to the speech 
detection unit 61 and a separate speech preprocessor may 
be used to parameterise the input speech for use by the 
speech recognition unit 25- However, such separate 
parameterisation of the input speech is not preferred 
because of the additional processing overhead involved. 

In the first embodiment described above, the AR filter 
coefficients output by the statistical analysis unit 21 
were converted into cepstral coefficients since the 
speech recognition unit used in the first embodiment was 



a cepstral based system. As those skilled in the art will 
appreciate, if the speech recognition system is designed 
to work with other spectral coefficients, then the 
coefficient converter 2 3 may be arranged to convert the 
AR filter coefficients into the appropriate spectral 
parameters. Alternatively still, if the speech 

recognition system is designed to operate with AR 
coefficients, then the coefficient converter 2 3 is 
unnecessary. 

In the above embodiments, Gaussian and Inverse Gamma 
distributions were used to model the various prior 
probability density functions of equation (19). As those 
skilled in the art of statistical analysis will 
appreciate, the reason these distributions were chosen is 
that they are conjugate to one another. This means that 
each of the conditional probability density functions 
which are used in the Gibbs sampler will also either be 
Gaussian or Inverse Gamma. This therefore simplifies the 
task of drawing samples from the conditional probability 
densities. However, this is not essential. The noise 
probability density functions could be modelled by 
Laplacian or student-t distributions rather than Gaussian 
distributions. Similarly, the probability density 
functions for the variances may be modelled by a 
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distribution other than the Inverse Gamma distribution. 
For excimple, they can be modelled by a Rayleigh 
distribution or some other distribution which is always 
positive- However, the use of probability density 
5 functions that are not conjugate will result in increased 

complexity in drawing samples from the conditional 
densities by the Gibbs sampler. 

Additionally, whilst the Gibbs sampler was used to draw 
10 samples from the probability density function given in 

equation (19), other seunpling algorithms could be used. 
For example the Metropolis-Hastings algorithm (which is 
reviewed together with other techniques in a paper 
entitled "Probabilistic inference using Markov chain 
15 Monte Carlo methods" by R. Neal, Technical Report CRG- 

TR-93-1, Department of Computer Science, University of 
Toronto, 1993) may be used to seunple this probability 
density - 

2 0 In the above embodiment, a Simulation Smoother was used 

to generate estimates for the raw speech samples. This 
Simulation Smoother included a Kalman filter stage and a 
smoothing filter stage in order to generate the estimates 
of the raw speech samples- In an alternative embodiment, 

25 the smoothing filter stage may be omitted, since the 
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Kalman filter stage generates estimates of the raw speech 
(see equation (33)). However, these raw speech samples 
were ignored, since the speech samples generated by the 
smoothing filter are considered to be more accurate and 
robust. This is because the Kalman filter essentially 
generates a point estimate of the speech samples from the 
joint probability density function P( s(n) | a^Jc^Oe^) , 
whereas the Simulation Smoother draws a sample from this 
probability density function. 

In the above embodiment, a Simulation Smoother was used 
in order to generate estimates of the raw speech samples. 
It is possible to avoid having to estimate the raw speech 
samples by treating them as "nuisance parameters" and 
integrating them out of equation (19). However, this is 
not preferred, since the resulting integral will have a 
much more complex form than the Gaussian and Inverse 
Gamma mixture defined in equation (19). This in turn will 
result in more complex conditional probabilities 
corresponding to equations (20) to (30). In a similar 
way, the other nuisance parameters (such as the 
coefficient variances or any of the Inverse Gamma, alpha 
and beta parameters) may be integrated out as well. 
However, again this is not preferred, since it increases 
the complexity of the density function to be sampled 
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using the Gibbs saunpler. The technique of integrating 
out nuisance parameters is well known in the field of 
statistical analysis and will not be described further 
here . 

5 

In the above embodiment, the data analysis unit analysed 
the samples drawn by the Gibbs sampler by determining a 
histogram for each of the model parameters and then 
determining the value of the model parameter using a 

10 weighted average of the samples drawn by the Gibbs 

sampler with the weighting being dependent upon the 
number of samples in the corresponding bin. In an 
alterative embodiment, the value of the model parameter 
may be determined from the histogram as being the value 

15 of the model parameter having the highest count. 

Alternatively, a predetermined curve (such as a bell 
curve) could be fitted to the histogram in order to 
identify the maximum which best fits the histogram. 

20 In the above embodiment, the statistical analysis unit 

modelled the underlying speech production process with a 
separate speech source model (AR filter) and a channel 
model. Whilst this is the preferred model structure, the 
underlying speech production process may be modelled 

25 without the channel model. In this case, there is no 



need to estimate the values of the raw speech samples 
using a Kalman filter or the like, although this can 
still be done. However, such a model of the underlying 
speech production process is not preferred, since the 
speech model will inevitably represent aspects of the 
channel as well as the speech. Further, although the 
statistical analysis unit described above ran a model 
order selection routine in order to allow the model 
orders of the AR filter model and the channel model to 
vary, this is not essential. In particular, the model 
order of the AR filter model and the channel model may be 
fixed in advance, although this is not preferred since it 
will inevitably introduce errors into the representation. 

In the above embodiments, the speech that was processed 
was received from a user via a microphone. As those 
skilled in the art will appreciate, the speech may be 
received from a telephone line or may have been stored on 
a recording medium. In this case, the channel model will 
compensate for this so that the AR filter coefficients 
representative of the actual speech that has been spoken 
should not be significantly affected. 

In the above embodiments, during the running of the model 
order selection routine, a new model order was proposed 



by drawing a random variable from a predetermined 
Laplacian distribution function. As those skilled in 
the art will appreciate, other techniques may be used. 
For example the new model order may be proposed in a 
deterministic way (ie under predetermined rules), 
provided that the model order space is sufficiently 
sampled. 
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CLAIMS 

1- A speech processing apparatus comprising: 

means for receiving a set of signal values 
5 representative of a speech signal generated by a speech 

source as distorted by a transmission channel between the 
speech source and the receiving means; 

a memory for storing data defining a predetermined 
function derived from a predetermined signal model which 
10 includes a first part having first parameters which 

models said source and a second part having second 
parameters which models said channel and said function 
being in terms of said first and second parameters and in 
terms of a set of signal values; 
15 means for applying said set of received signal 

values to said stored function; and 

means for processing said function with those values 
applied to obtain values of said first parameters that 
are representative of said speech generated by said 
20 speech source before it was distorted by said 

transmission channel. 



25 



2. A speech processing apparatus comprising: 

means for receiving a set of signal values 
representative of a speech signal generated by a speech 



63 

source as modified by a transmission channel between the 
speech source and the receiving means; 

means for storing data defining a predetermined 
function derived from a predetermined signal model which 
5 includes a first part having first parameters which 

models said source and a second part having second 
parameters which models said channel, said function being 
in terms of said first and second parameters and 
generating, for a given set of signal values, a 

10 probability density function which defines, for a given 

set of first parcuneters and second parameters, the 
probability that the predetermined signal model has those 
parameter values, given that the signal model is assumed 
to have generated the received set of signal values; 

15 means for applying said set of received signal 

values to said function; 

means for processing said function with those values 
applied to derive scimples of at least said first 
parameters from said probability density function; and 

20 means for analysing at least some of said derived 

samples to determine values of said first pareimeters that 
are representative of said speech signal generated by 
said source before it was modified by said transmission 
channel . . 



25 



3. A speech processing apparatus comprising: 

means for receiving a set of signal values 
representative of a speech signal generated by a speech 
source as distorted by a transmission channel between the 
speech source and the receiving means; 

a memory for storing data defining a predetermined 
function derived from a predetermined signal model which 
includes a first part having first parameters which 
models said source and a second part having second 
parameters which models said channel, said function being 
in terms of said first and second parameters and being in 
terms of a set of raw speech signal values representative 
of speech generated by said source before being distorted 
by said transmission channel; 

means for processing said received set of signal 
values with initial estimates of said first and second 
parameters, to generate an estimate of the raw speech 
signal values corresponding to the received set of signal 
values ; 

means for applying said set of received signal 
values and the estimated set of raw speech signal values 
to said function; 

means for processing said function with those values 
applied to obtain values of said first parameters that 
are representative of said speech signal generated by 
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said speech source before it was distorted by said 
transmission channel. 

4- A speech processing apparatus comprising: 
5 means for receiving a set of signal values 

representative of a speech signal generated by a speech 
source as modified by a transmission channel between the 
signal source and the receiving means; 

a memory for storing data defining a predetermined 

10 function derived from a predetermined signal model which 

includes a first part having first parameters which models 
said source and a second part having second pareuneters 
which models said channel, said function being in terms 
of said first and second parameters and being in terms 

15 of a set of raw speech signal values representative of a 

speech signal generated by said source before being 
modified by said transmission channel, said function 
generating, for a given set of signal values, . a 
probability density function which defines, for a given 

20 set of first parameters, second parameters and raw speech 

signal values, the probability that the predetermined 
signal model has those parsimeter values and generates said 
raw speech signal values, given that the signal model is 
assumed to have generated the received set of signal 

25 values; 
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means for processing said received set of signal 
values with initial estimates of said first and second 
parameters, to generate an estimate of said set of raw 
speech signal values corresponding to the received set of 
5 signal values; 

means for applying said set of received signal values 
and the estimated set of raw speech signal values to said 
function; 

means for processing said function with those values 
10 applied to derive samples of at least said first 

parameters from said probability density function; and 

means for analysing at least some of said derived 
samples to determine values of said first pareuneters that 
are representative of said speech signal generated by said 
15 speech source before it was modified by said transmission 

channel . 



5 . An apparatus for determining sets of paraimeter values 
representative of an input speech signal, the apparatus 
20 comprising: . . 

means for receiving a plurality of speech signal 
values representative of an input speech signal; 

means for dividing the plurality of speech signal 
values into a succession of groups of speech signal 
25 values; and 



means for determining a set of parameter values 
representative of the speech signal values in each group; 

wherein said determining means comprises means for 
varying the number of parameter values used to represent 
the speech signal values in each group. 

6. An apparatus for determining sets of parameter values 
representative of an input speech signal, the apparatus 
comprising: 

means for receiving a plurality of speech signal- 
values representative of an input speech signal; 

means for dividing the plurality of speech signal 
values into a succession of groups of speech signal,, 
values; and 

means for processing the speech signal values in each 
group to determine a set of parameter values 
representative of the speech signal values in the group; 
wherein said processing means comprises: 
a memory for storing data defining a predetermined 
function which gives, for a set of speech signal values 
of a group, a probability density for parameters of a 
predetermined signal model which is assumed to have 
generated the speech signal values in the group, the 
probability density defining, for a given set of parameter 
values, the probability that the predetermined signal 
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model has those parameter values, given that the model is 
assumed to have generated the speech signal values in the 
group; 

means for applying the set of speech signal values 
of a current group to said stored function to give the 
probability density for said model parameters for the 
current group; 

means for processing said function to derive samples 
of parameter values from said probability density for the 
current, group; 

means for evaluating said probability density for the 
current group using one or more of said derived samples 
of parameter values for different numbers of parameter 
values to determine respective probabilities that the 
predetermined signal model has those parameter values; and 

means for processing at least some of said derived 
Scunples of parameter values and said evaluated 
probabilities to determine model parameters that are 
representative of the set of signal values in the current 
group. 

7. An audio comparison apparatus comprising: 

a memory for storing a predetermined function which 
gives, .for a given set of audio signal values, a 
probability density for parameters of a predetermined 






5 



10 



15 



20 



audio model which is assumed to have generated the set of 
audio signal values, the probability density defining, for 
a given set of model parameter values, the probability 
that the predetermined audio model has those parameter 
values, given that the model is assumed to have generated 
the set of audio signal values; 

means for receiving a set of audio signal values 
representative of an input audio signal; 

means for applying the set of received audio signal 
values to said stored function to give the probability 
density for said model parameters for the set of received 
audio signal values; 

means for processing said function, with said set.. of 
received audio signal values applied, to derive samples 
of parameter values from said probability density; 

means for analysing at least some of said derived 
Scimples of parcimeter values to determine parameter values 
that are representative of the set of received audio 
signal values; and 

means for comparing said determined parameter values 
with pre-stored parameter values to generate a comparison 



result . 



8. 



A speech recognition system comprising an audio 



comparison apparatus according to claim 7 . 
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9- A speaker verification system comprising an audio 
comparison apparatus according to claim 7 . 

10- An acoustic classification system comprising an audio 
5 comparison apparatus according to claim 7 . 

11. An audio transmission system comprising: 

a memory for storing a predetermined function which 
gives, for a given set of audio signal values, a 
10 probability density for parameters of a predetermined 

audio model which is assumed to have generated the set of 
audio signal values, the probability density defining, for 
a given set of model parameter values, the probability 
that the predetermined audio model has those parameter 
15 values, given that the model is assumed to have generated 

the set of audio signal values; 

means for receiving a set of audio signal values 
representative of an input audio signal; 

means for applying the set of received audio signal 
2 0 values to said . stored function to give the probability 

density for said model parameters for the set of received 
audio signal values; 

means for processing said function with said set of 
received audio signal values applied, to derive scimples 
25 of parameter values from said probability density; 
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means for analysing at least some of said derived 
samples of parameter values to determine parameter values 
that are representative of the set of received audio 
signal values; 

5 means for transmitting said determined pareimeter 

values; 

means for receiving said transmitted parameter 
values ; and 

means for processing the received parameter values 
10 to generate an output signal in dependence thereon- ^ 

12. An apparatus for determining a quality measure 
indicative of the quality of an audio signal, the 
apparatus comprising: 

15 a memory for storing a predetermined function which 

gives, for a given set of audio signal values, a 
probability density for parameters of a predetermined 
audio model which is assumed to have generated the set of 
audio signal values, the probability density defining, for 

2 0 a given set of model parameter values, the probability 

that the predetermined audio model has those parameter 
values, given that the model is assumed to have generated 
the set of audio signal values; 

means for receiving a set of audio signal values 

2 5 representative of an input audio signal; 



means for applying the set of received audio signal 
values to said stored function to give the probability 
density for said model parameters for the set of received 
audio signal values; 

means for processing said function with said set of 
received audio signal values applied, to derive samples 
of parameter values from said probability density; 

means for analysing at least some of said derived 
samples of parameter values to determine a measure of the 
variance of said samples of parameter values; and 

means for outputting a signal indicative of the 
quality of the received audio signal values independence 
upon said determined variance measure. 

13 . An audio processing system substantially as 
hereinbefore described with reference to or as shown in 
the accompanying drawings. 

14. An audio processing method substantially as herein 
above described with reference to or as shown in the 
accompanying drawings . 

15 . A computer readable medium storing computer 
executable process steps to configure a programmable 
computer apparatus to be configured as an apparatus 



according to any of claims 1 to 12 . 



16. A processor implement able process steps for 

configuring a programmable computing device to be 
configured as an apparatus according to any of claims 1 
to 12. 



ABSTRACT 



SPEECH PROCESSING SYSTEM 



A speech processing system is provided which is operable 
to receive sets of signal values representative of a 
speech signal generated by a speech source as distorted 
by a transmission channel between the speech source and 
the speech processing system. The system stores data 
defining a predetermined function derived from a signal 
model which models both the speech source and the channel 
and defining a probability density function which gives, 
for a given set of model parameters, the probability that 
the signal model has those model pareoneters given that 
the signal model is assumed to have generated the 
received set of signal values. The system applies a 
current set of received signal values to the stored 
probability density function and then draws samples from 
it using a Gibbs sampler. The system then analyses the 
samples to determine a set of parameter values 
representative of the speech signal before it was 
distorted by the channel. 




FIG. 1 



o 

LU 
LU 


DELS 


Q- 


O 


CO 





CM 




CO 
CM" 



CM' 



CJl 



I— 


cc 


UJ 


LU 


CD 


1— 
DC 
LU 


LL. 
U- 
UJ 


> 


O 


o 


O 





<C CO 

CO 




3/0 

FIG. 4 



0 



START 



S1' 



0 



PROPOSE A NEW 
MODEL ORDER (K2) 



S3 



CALCULATE: 



MO = Min 



P(a<i:K3>.'<2l ) 

P(g<l ;K,>>'<i I ) 



.1 



YES 




NO 



SET MODEL ORDER TO 
K2 AND INCREMENT A 
COUNT ASSOCIATED 
WITH VALUE OF K2 



T 

S7 



LEAVE MODEL ORDER AT 
CURRENT VALUE (Ki) AND 
INCREMENT A COUNT 
ASSOCIATED WITH 
VALUE OF Ki 



S9 



A-/13 

(^START^ 



F\G. 5 



S29 



i 



INCREMENT t 



S37 



DECREMENT t 



INITIALISE TIME 
VARIABLE t 



S21 



COMPUTE KALMAN FILTER 
EQUATIONS FOR SAMPLES 
AT CURRENT TIME (t) 



S23 



STORE KALMAN FILTER 
VARIABLES COMPUTED 
FOR CURRENT TIME (t) 



S25 




S27 



COMPUTE SMOOTHING FILTER 
EQUATIONS USING KALMAN 
FILTER VARIABLES COMPUTED 
FOR CURRENT TIME(t) 



331 





CM 1 



CO 



<£* CO 



>- 
o 







•< 


:^ 


DC 




> 


cc 







L 




oc o 



Q UJ 

o ^ 









oc 




CD 




lOOTHI 


LMAN 
LTER 


OTHIN 
LTER 






<c u: 


O CE 




CO 












CO 




o 












( 




1 




ZD 




1 




CO 


CO 







00 



FIG. 7 



y,{i) — ' y|(N) 



a9, 


h9, 




r9 






C'e ■ 




PI' 










e 







forg = 0 to M 



w{t). d(t). L(t) 
fort = 1 to N 



s^(t) and e^(t) 
for t = n-v-N+1 to n 
where v = max (k.r) 



k MODEL ORDER COUNTS 
r MODEL ORDER COUNTS 



1 



(^START^ 



7/13 

FIG. 8 



RETRIEVE NEXT FRAME OF SPEECH 
SAMPLES FROM THE BUFFER 



S41 



SET INITIAL ESTIMATES FOR 
CHANNEL MODEL.SPEECH MODEL 
AND SPEECH AND NOISE STATISTICS 

i 



S43 



USE SIMULATION SMOOTHER 
TO PROVIDE ESTIMATE OF 
RAW SPEECH SAMPLES 



S45 



USE ESTIMATED SPEECH SAMPLES 
AND GIBBS SAMPLER TO UPDATE 
CHANNEL MODEL AND SPEECH 
MODEL AND NOISE STATISTICS 



S49 



S47 

UPDATE' 
MODEL ORDERS. 
? 

,YES 



NO. 



1JPDATE 
RAW SPEECH 
.SAMPLES^ VES 

" S55 



PERFORM MODEL ORDER 
SELECTION USING REVERSIBLE 

JUMP MARKOV CHAIN 
MONTE CARLO C OMPUTATION 

T" 

S51 



S53 

ANOTHER" 
"GIBBS iteration: 
? 

NO 



YES 



S57 
i 



DETERMINE BEST 
MODEL ORDERS 







OUTPUT CHANNEL 
MODEL.SPEECH 
MODEL AND 
SPEECH AND NOISE 
STATISTICS FOR 
CURRENT FRAME 

\ 



S59 




YES 



S61 



FIG. 9a 



COUNT 



1 r 

2 3 4 5 



T r 
6 7 8 



1 r 



AR FILTER MODEL ORDER (k) 



COUNT 



0 
0.26 



FIG. 9b 



Inn, 



0.28 



0.3 



0.32 



0.34 



0.36 



PROCESS NOISE VARIANCE (o^ ) 



COUNT 



FIG. 9c 



0' 

-0.7 -0.65 -0.6 -0.55 -0.5 -0.45 -0.4 



o 

C/5 S 







o 






l-j 




ID 


CD 


CO 


ICO 


LU 

or 


1 1 1 

DC 


A 



V 

CM 



o 

Si 

UJ ^ 

LU CD 
Q_ o 
CO o 

LU 

cc 



CO 



□C CD 

o p: != 

LU ^ 

l^^ ^ § 

Q_ I— —J 

CO LU 
O 




CD 



CM" 



CD 



■< CO 

li 

CO 



CQ 




CD 



CD 



•< CO 
CO 



c 

L 


c 


LL. 
CQ 


y 




A/D 


y 




a: 

LU 


r 
L 


L. 



nli3 



u- tn 

CO ^ 
CO 

o 



1 







ACOUSTIC 
LASSIFICATION 
SYSTEM 


ISTIC 

IcLS 


■< 










o 



cvj 
CD 




nli3 




eg- 



os 



r 

LU 


DC 




LU 




O 




CJ) 






CD 


LU 



<C CO 



0 

u 

t 


C 
u 


■ ' 

Lt_ 
ZD 
CQ 


i 




A/D 


/ 




CC 
LU 


r 

ij 


I- 







S£2 










o 


LU 




LU 
LU 
CU 


HiN 


UNI 


C/D 








CO 





CD <^ 
O " ' 



CD 
LU 

DC 



CC 



00 




/J//3 



CO 
<C 
CO 



!< =! 



V 

CD 



CD 



UJ o 

^ 5 t= 



^ cl Q 
cj) t: CO 

LU _J CO 
LU -=C LU 

CI- CO 
CO o ^ 

^ 



<C CO 

CO 

CO 




DC 



CD 



DC 





FITZPATRICK, CELLA, HARPER, & SCINTO 

30 Rockefeller Plaza 
New York, NY 10112-3801 



